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In two dimensions, quenched disorder always rounds transitions involving the breaking of spatial 
symmetries so, in practice, it can often be difficult to infer what form the symmetry breaking would 
take in the "ideal," zero disorder limit. We discuss methods of data analysis which can be useful 
for making such inferences, and apply them to the problem of determining whether the preferred 
order in the cuprates is "stripes" or "checkerboards." In many cases we show that the experiments 
clearly indicate stripe order, while in others (where the observed correlation length is short), the 
answer is presently uncertain. 



I. INTRODUCTION 



Charge ordered states are common in strongly cor- 
related materials, including especially the cuprate high 
temperature superconductors. Identifying where such 
phases occur in the phase diagram, and where they oc- 
cur as significant fluctuating orders is a critical step in 
understanding what role they play in the physics, more 
generally. Since "charge ordered" refers to states which 
spontaneously break the spatial symmetries of the host 
crystal, identifying them would seem to be straightfor- 
ward. However, two real- world issues make this less sim- 
ple than it would seem. In the first place, quenched dis- 
order (alas, an unavoidable presence in real materials), 
in all but a very few special circumstances, rounds the 
transition and spoils any sharp distinction between the 
symmetric and broken symmetry states. Moreover, the 
charge modulations involved tend to be rather small in 
magnitude, and so difficult to detect directly in the obvi- 
ous experiments, such as X-ray scattering. 

In a previous paper ji three of us addressed at some 
length the issue of how the presence or absence of charge 
order or incipient charge order can best be established 
in experiment. In the present paper we focus on a re- 
lated issue: in a system in which charge order is believed 
to exist, how can the precise character of the charge 
order best be established? This is particularly timely 
given the spectacular developments in scanning tunneling 
microscopy (STM) which produces extremely evocative 
atomic scale "pictures" of the local electronic structure 
- the question is how to extract unambiguous conclu- 
sions from the cornucopick2iiiiSiSiL2iS of data. We take 
as a representative example, the issue of whether the 
charge order that is widely observed in the cuprates is 
"stripes" (which in addition to breaking the translation 
symmetry, breaks various mirror and discrete rotation 
symmetries of the crystal) or "checkerboards" (an order 
which preserves the point-group symmetries of the crys- 
tal). To address this issue, we generate simulated STM 
data and then test the utility of various measures we have 
developed for discriminating different types of order by 



applying them to this simulated data. Where the corre- 
lation length for the charge order is long, definitive con- 
clusions can be drawn relatively simply - consequently, 
it is possible to conclude that the preferred charge or- 
der in the 214 (La2Cu04) family of materials is stripes 
and not checkerboards However, where the correlation 
length is short (disorder effects are strong), it turns out 
(unsurprisingly) to be very difficult to develop any fool- 
proof way to tell whether the observed short-range order 
comes from pinned stripes or pinned checkerboards - for 
example, the image in Fig. ^ (right panel) corresponds 
to disorder- pinned stripes, despite the fact that, to the 
eye, the pattern is more suggestive of checkerboard order 
(with the latter seen in Fig. [3 (right panel)). 

In Section^ we give precise meaning in terms of bro- 
ken symmetries to various colloquially used descriptive 
terms such as "stripes," "checkerboards," "commensu- 
rate," "incommensurate," "diagonal," "vertical," "bond- 
centered," and "site-centered." In Section lTlll we write an 
explicit Landau-Ginzburg (LG) free energy functional for 
stripe and checkerboard orders, including the interactions 
between the charge order and impurities. In Section llVI 
we generate simulated STM data by minimizing the LG 
free energy in the presence of disorder. (See Figs. ^ and 
121) The idea is to develop strategies for solving the in- 
verse problem: Given the simulated data, how do we 
determine whether the "ideal" system, in the absence of 
disorder, would be stripe or checkerboard ordered, and 
indeed, whether it would be ordered at all or merely in 
a fiuctuating phase with a large CDW susceptibility re- 
fiecting the proximity of an ordered state. In Section Hvl 
we define several quantitative indicators of orientational 
order that are useful in this regard, but unless the cor- 
relation length is well in excess of the CDW period, no 
strategy we have found allows confident conclusions. In 
Section |Vl we show that the response of the CDW or- 
der parameter to various small symmetry breaking fields, 
such as a small orthorhombic distortion of the host crys- 
tal, can be used to distinguish different forms of charge 
order. In Section IVII we apply our quantitative indica- 
tors to a sample of STM data in Bi2Sr2CaCu208+5 and 
discuss the results. In Section IVIII we conclude with a 
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few general observations. 



II. GENERAL CONSIDERATIONS 

Stripes are a form of unidirectional charge order (see 
Fig. n (left panel), characterized by modulations of the 
charge density at a single ordering vector, Q, and its 
harmonics, Q„ ~ nQ with n = an integer. In a crys- 
tal, we can distinguish different stripe states not only 
by the magnitude of Q, but also by whether the order 
is commensurate (when |Q|a = 27r(m/n) where a is a 
lattice constant and n is the order of the commensurabil- 
ity) or incommensurate with the underlying crystal, and 
on the basis of whether Q lies along a symmetry axis 
or not. In the cuprates, stripes that lie along or nearly 
along the Cu-0 bond direction are called "vertical" and 
those at roughly 45° to this axis are called "diagonal." 
In the case of commensurate order, stripes can also be 
classified by differing patterns of point-group symmetry 
breaking - for instance, the precise meaning of the often 
made distinction between so-called "bond-centered" and 
"site-centered" stripes is that they each leave different 
reflection planes of the underlying crystal unbroken. Fur- 
thermore, it has been argued that bond and site-centered 
stripes may be found in the same material, and even 
may coexist at the same temperatureji^ The distinction 
between bond and site-centered does not exist for incom- 
mensurate stripes. If the stripes are commensurate, then 
Q must lie along a symmetry direction, while if the CDW 
is incommensurate, it sometimes will not. 

Checkerboards are a form of charge order (see Fig. [21 
(left panel) that is characterized by bi-directional charge 
density modulations, with a pair of ordering vectors, Qi 
and Q2 (where typically |Qi| = |Q2|-) Checkerboard or- 
der generally preserves the point group symmetry of the 
underlying crystal if both ordering vectors lie along the 
crystal axes. In the case in which they do not, the or- 
der is rhombohedral checkerboard and the point group 
symmetry is not preserved. As with stripe order, the 
wave vectors can be incommensurate or commensurate, 
and in the latter case Qj a = 27r (m/rt, m'/n'). Com- 
mensurate order, as with stripes, can be site-centered or 
bond-centered. 



III. LANDAU-GINZBURG EFFECTIVE 
HAMILTONIAN 

To begin with, we will consider an idealized two di- 
mensional model in which we ignore the coupling be- 
tween layers and take the underlying crystal to have the 
symmetries of a square lattice. We further assume that 
in the possible ordered states, the CDW ordering vector 
lies along one of a pair of the orthogonal symmetry di- 
rections, which we will call "x" and "y". We can thus 
describe the density variations in terms of two complex 
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TABLE L Phases of the Landau-Ginzburg model, in the ab- 
sence of disorder. 

scalar order parameters, 

p{r) = p+yi{r)e'^-'' + ip2{r)e'^^y + C.C.] (3.1) 

(For simplicity, we will take Qx = Qy throughout.) Note, 
the "density, " in this case, can be taken to be any scalar 
quantity, for instance the local density of states, and need 
not mean, exclusively, the charge density. 

To quartic order in these fields and lowest order in 
derivatives, and assuming that commensurability effects 
can be neglected, the most general Landau-Ginzburg ef- 
fective Hamiltonian density consistent with symmetry 
has been written down by several authors 

(3.2) 

The sign of a determines whether one is in the broken 
symmetry phase (a < 0) or the symmetric phase {a > 0), 
and in the broken symmetry phase, 7 determines whether 
the preferred order is stripes (7 > 0) or checkerboards 
(7 < 0). Note that for stability, it is necessary that 
7 > ~u and u > 0; if these conditions are violated, one 
needs to include higher order terms in H^ff- Without 
loss of generality, we can rescale distance so that kl = 1 
and the order parameter magnitude such that u = 1. For 
simplicity, in the present paper, we will also set kl = kt, 
although the more general situations can be treated with- 
out difficulty. The phases of this model in the absence of 
disorder are summarized in Table 

Imperfections of the host crystal enter the problem as 
a quenched potential, U{r): 

Hd,s - C/(r)p(r) (3.3) 

To be explicit, we will take a model of the disorder 
potential in which there is a concentration of impuri- 
ties per unit area, S/a^ where a is the "range" of the 
impurity potential and Uq is the impurity strength, so 
Ui'r) = X)i UqQIo^ ^ (r — fj)^]; where the sum is over 
the (randomly distributed) impurity sites, and Q is 
the Heaviside function. We have arbitrarily taken a to 
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be 1/4 the period, A, of the CDW, i.e. A = 27r/|Q| and 
a — A/4. (This choice is motivated by the fact that, in 
many cases, the observed charge order has a period A w 4 
lattice constants.) 



IV. ANALYSIS OF THE SIMULATED DATA 

In this section we will show how these ideas can be used 
to interpret STM images in terms of local stripe order. 
In Ref.|ll| it was shown that local spectral properties of 
the electron Green function of a correlated electron sys- 
tem, integrated over an energy range over a window in 
the physically relevant low energy regime, can be used as 
a measure of the local order. This is so even in cases in 
which the system is in a phase without long range order 
but close enough to a quantum phase transition ( "fluctu- 
ating order") that local defects can induce local patches 
of static order. From this point of view any experimen- 
tally accessible probe with the correct symmetry can be 
used to construct an image of the local order state. In 
applying the following method to real experimental data, 
one must take as a working assumption that the image 
obtained is representative of some underlying order, be it 
long-ranged or incipient. This analysis, of course, would 
not make sense if the data is not, at least in substan- 
tial part, dominated by the correlations implied by the 
existence of an order parameter. 

We generate simulated data as follows: For a given 
randomly chosen configuration of impurity sites, we min- 
imize -ffe// + Hj^is with respect to if. This is done nu- 
merically using Newton's method. The order parameter 
texture is then used to compute the resulting density map 
according to Eq. 13.11 This we then treat as if it were the 
result of a local imaging experiment, such as an STM 
experiment. 

Even weak disorder has a profound effect on the re- 
sults. For a < 0, collective pinning causes the broken 
symmetry state to break into domains with a charac- 
teristic size which diverges exponentially as fJo ~> (In 
three dimensions, the ordered state survives as long as 
the disorder is less than a critical value.) Examples of 
this are shown in Fig. ^ and Fig. [21 where data with 
a given configuration of impurities with concentration 
d = 0.1 are shown for various strengths of the poten- 
tial, Uq. For a checkerboard phase (7 < 0), the domain 
structure is rather subtle, involving shifts of the phase of 
the density wave as a function of position as can be seen 
in Figs. |2Ileft and center panels). In the stripe phase, in 
addition to phase disorder, there is a disordering of the 
orientational ("electron nematic") order, resulting in a 
more visually dramatic breakup into regions of vertical 
and horizontal stripes, as can be seen in Fig. ^ (center 
panel). 

The effect of quenched disorder in the symmetric phase 
{a > 0) is somewhat different. In a sense, the effect of the 
disorder is to pin the fluctuating order of the proximate 
ordered phase. However, here, whether the disorder is 



weak or strong, it is nearly impossible to distinguish fluc- 
tuating stripes from checkerboards. Fig. ^ (right panel) 
and Fig. 13 (right panel) illustrate this phenomenon. This 
is easily understood in the weak disorder limit, where 

ifijir) = J dr'xo(r - r')e-^'J"-'c/(r') -I- ©(C/^) (4.1) 

where the susceptibility, 

Xoir) = KoiV^r), (4.2) 

is expressed in terms of the Kq Bessel function and is 
independent of 7. Near criticality (1 ^ a > 0), the 
susceptibility is very long ranged, so a significant degree 
of local order can be pinned by even a rather weak im- 
purity potential. However, only the higher order terms 
contain any information at all about the sign of 7, and 
by the time they are important, the disorder is probably 
already so strong that it blurs the distinction between 
the two states, anyway. 



A. Diagnostic Filters 

Now, our task is to answer the question: Given a set 
of simulated data, what quantitative criteria best allows 
us to infer the form of the relevant correlations in the 
absence of disorder? For sufficiently weak disorder, these 
criteria are, at best, just a way of quantifying a con- 
clusion that is already apparent from a visual analysis 
of the data. Where disorder is of moderate strength, 
such criteria may permit us to reach conclusions that are 
somewhat less prejudiced by our preconceived notions. 
Of course, when the disorder is sufficiently strong that 
the density-wave correlation length is comparable to the 
CDW period, it is unlikely that any method of analysis 
can yield a reliable answer to this question. 

Firstly, to eliminate the rapid spatial oscillations, we 
define two scalar fields (which we will consider to be the 
two components of a vector field, A(r)) corresponding 
to the components of the density which oscillate, respec- 
tively, with wave vectors near Qx and Qy: 

A,(r)=y'dr'F,(r-r')p(r') (4.3) 

where we take Fj to be the coherent state with spatial 
extent equal to the CDW period: 

^j(r) = ^ exp[z Q,r, - rV(2^A2)] (4.4) 

(no summation over j.) 

In terms of A we construct three quantities which can 
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FIG. 1: (color online) left panel: Highly stripe-ordered system, with weak impurities, Uo = 0.1, 6 = 0.1. Here 7 = 1, a = —0.05. 
[Scale is arbitrary.] center panel: Otherwise identical to the first system (including the spatial distribution and concentration 
of impurities), but the strength each impurities has increased to Uo ~ 0.75. right panel: Identical to the left panel, except 
a = +0.05. Much of the underlying charge pattern remains, even to positive a, where in the absence of impurities, the system 
would be homogeneous. All graphs are approximately 20 CDW wavelengths in width. 
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FIG. 2: (color online) The parameters entering the effective Hamiltonian and the impurity realizations are identical here to the 
panels of Fig. Q with the exception of the symmetry breaking term, 7, which is now —0.95. (In the center panel, because the 
checkerboard state is more stable than the analogous stripe state, we have taken Uo — 1.5.) Unlike the stripe ordered system, 
the checkerboard system does not break into domains, but rather develops pair wise dislocations. In |21 (center panel), three 
pairs of such dislocations are visible. Note the similarity between the right panel of each set of Fig^and Fig|5| the sign of 7 
has little effect for q > 0. 



be used in interpreting data: 
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(4.6) 



The quantities called ^ have units of length and rjorient is 



dimensionless. All of these quantities are invariant under 
a change of units, p{r) — > Ap(r); this is important since 
in many experiments, including STM, the absolute scale 
of the density oscillations is difficult to determine because 
of the presence of unknown matrix elements. 

^CDW has the interpretation of a CDW correlation 
length. In the absence of quenched disorder, and for 
a < 0, £,CDW ~ where L is the linear dimension of 
the sample. In the presence of disorder, ^cdw is an av- 
erage measure of the domain size. For a > and weak 



(4.7) but non- vanishing disorder, £,cdw 



as can be 



seen from a scaling analysis of Eq. 14.11 The evolution of 
^CDW as a function of a is shown in Fig.|31for a system of 
size L = 20A, for various strengths of the disorder and for 
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FIG. 3: (color online) ^cdw vs. a (measured in units of 
the CDW wavelength), top panel) 7 = 1, bottom panel) 7 = 
—0.95. In a perfectly clean system, (,cdw vanishes for a > 0, 
whereas with even a little disorder, charge order is induced. 
For Uo > 1 and a < 0, disorder affects ^cdw more strongly 
in the stripe system. For a > 0, there is little distinction in 
either the sign of 7 or the strength of Uq. All quantities in 
Figs. I^m and 1^ are computed for systems of size 20 A x 20 A 
and averaged over 50 or more realizations of the disorder. 



stripes (Fig.|31(top panel) and checkerboards (Fig. O (bot- 
tom panel).) Ccdw is generally a decreasing function of 
increasing disorder, although for a > there is a range 
in which it exhibits the opposite behavior. For fixed, 
non-zero disorder, we see that a large value oi ^cdw > 4 
almost inevitably means that a < 0, i.e. that the den- 
sity patterns are related to a domain structure of what 
would otherwise have been a fully ordered state. How- 
ever, smaller values of S,cdw can either come from weak 
pinning of CDW order which would otherwise be in a 
fluctuating phase, or a very small domain structure due 
to strong disorder. 

The CDW correlation length does not distinguish be- 
tween stripe and checkerboard patterns. However, for 
a < 0, the orientational amplitude rjorient is an effec- 
tive measure of stripiness. In the clean system with 
a < 0, rjorient approaches unity for 7 > and is zero for 
7 < 0. While quenched disorder somewhat rounds the 
sharp transition in rjorient at 7 = 0, it is clear from Fig.^ 



(top panel) that values of rjorient > 0.2 are clear indica- 
tors of stripe order, and rjorient < 0.2 implies checker- 
board. In the absence of disorder, rjorient is ill-defined 
for a > 0, and even for non-zero disorder, the behavior 
of rjorient is difficult to interpret in the fluctuating order 
regime, as is also clear from Fig. 0] The orientational 
correlation length, £,orient, gives similar information as 
rjorient, and suffers from the same shortcomings. 

One interesting possibility is that, for a weakly disor- 
dered stripe phase, one can imagine an orientational glass 
in which £,orient ^ £,CDW, i-S- the CDW order is phase 
disordered on relatively short distances, but the orien- 
tational order is preserved to much longer distances. In 
Fig.O we plot the ratio of ^orient / £.c dw for 7 = 1 (strong 
preference for stripes) as a function of a for various val- 
ues of the disorder. Clearly, we have not found dramatic 
evidence of such an orientational glass, although we have 
not carried out an exhaustive search. Nonetheless, for 
a < 0, this ratio is manifestly another good way to dis- 
tinguish stripe and checkerboard order. 

The bottom line: If £,cdw is a few periods or more, 
it is possible to conclude that a < 0, i.e. that in the 
absence of impurities there would be long-range CDW 
order. If ^cdw is shorter than this, then either the impu- 
rity potential is very strong (which should be detectable 
in other ways) or a ~ ^cdw positive. For intermediate 
values of ^cdw, all that can be inferred is that the sys- 
tem is near critical, \a\ <C 1. Given a substantial CcdWj 
it is possible to distinguish a pinned stripe phase from 
a pinned checkerboard phase for which rjorient is greater 
than or less than 0.2, respectively. 



V. EFFECT OF AN ORTHORHOMBIC 
DISTORTION 

An orthorhombic distortion breaks the C4 symmetry 
of the square lattice down to C2. There are two distinct 
ways this can occur - either the square lattice can be 
distorted to form rectangles, as shown in Fig.|SJi, in which 
case the "preferred" orthorhombic axis is either vertical 
or horizontal, or the squares can be distorted to form 
rhombi, as shown in Fig. in which case the preferred 
orthorhombic axis is diagonal. A general orthorhombic 
distortion is represented by a traceless symmetric tensor, 
Oab] an orthorhombic distortion corresponding to Fig.|SlD 
is represented by O = has while Fig.|^ is O = hai where 
h is the magnitude of the symmetry breaking and aj are 
the Pauli matrices. Then 
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-OabQaQb [|(^1P-|(P2|'] 
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(5.1) 



where ... is higher order terms. 

In case (a), the first term is non-zero, and hence dom- 
inant. For h positive, this enhances Lpi and suppresses 
if 2. In a stripe phase, this has the same effect as a mag- 
netic field in a ferromagnet - it chooses among the oth- 
erwise degenerate vertical and horizontal stripe ordered 
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FIG. 4: (color online) rjorient vs. 7: top panel) a < 0. In the 
ordered phase, rjorient is good indicator of the nature of the 
underlying order (i.e the sign of 7. At large Uo, the distinc- 
tion is lost, and the result approaches that of the symmetric 
phase (a > 0), shown in the bottom panel. We observe that 
the nearly uniform value of rjorient ~ 0.2 in the a > mea- 
surements intersects the (all of the) data in the (a = —0.03) 
graph at the 7 = axis. 
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FIG. 5: (color online) £,orient/^CDW vs. a: For a < 0, 
iorient / £,c Dw is a stroug ludicator of the sign of 7. For a > 
and either sign of 7, the disorder-averaged ratio is 1/2, largely 
independent of other parameters. 



FIG. 6: Orthorhombic symmetry breaking reduces a square 
lattice to a lower symmetry, (a) Rectangular lattice distortion 
(exaggerated). The preferred orthorhombic axis lies along an 
original lattice vector (i.e. along the lines connecting atomic 
sites.) (b) A rhombohedral distortion leaves the preferred 
orthorhombic axis diagonal to the original lattice vectors. 



states, so one is preferrediii For checkerboard order, it 
produces a distortion of the fully ordered state, so that 
the expectation value of ipi exceeds the expectation value 
of ip2- Moreover, it results in a spUt phase transition, 
so that as a function of decreasing temperature, rather 
than a single transition from a symmetric high temper- 
ature phase to a low temperature checkerboard phase, 
in the orthorhombic case there are two transitions, the 
first to a stripe ordered phase, and then at a temper- 
ature smaller by an amount proportional to /i, a transi- 
tion to a distorted checkerboard phase. The second term, 
proportional to is subdominant in this case, but still 
has a significant effect. For an incommensurate stripe 
phase, it results in a small shift in the ordering wave 
vector Q ^ Q = Q(l — gh/ k.l)- In an incommensu- 
rate checkerboard phase, it results in a relative shift of 
the two ordering vectors, Q ^ Q = Q(l — gh/KL) and 
Q' ^ Q' = Q'(l + gh/nL) one toward smaller and the 
other toward larger magnitude producing a rectangular 
checkerboard. In the case in which the order is commen- 
surate, it is locked to the lattice, and therefore the only 
shifts in ordering wave-vectors are proportional to the 
(usually miniscule) shifts of the lattice constant. 

In case (b), the first term vanishes, so the second term 
is dominant. For incommensurate order, this results in 
a small rotation of the ordering vector away from the 
crystalline symmetry axis. To first order in ft,, the new 
ordering vector is Q = |Q|(l,fc) with k = gh/Kx and, 
in the case of checkerboard order, the second ordering 
vector is Q' = |Q|(fc, 1). Again, in the commensurate 
case, the order remains locked to the lattice until the 
magnitude of the orthorhombicity exceeds a finite critical 
magnitude. 

To summarize, the response of charge order to small 
amounts of orthorhombicity can be qualitatively differ- 
ent depending on whether the order is commensurate or 
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incommensurate and checkerboard or striped. 



1. More complex patterns of symmetry breaking 

It is useful to point out that with complex crystal 
structures, the application of the above ideas requires 
some care. For example, there are some cuprate su- 
perconductors which exhibit a so called Low Tempera- 
ture Tetragonal (LTT) phase. This phase has an effec- 
tive orthorhombic distortion of each copper oxide plane, 
but has two planes per unit cell and a four-fold twist 
axis which is responsible for the fact that it is classi- 
fied as tetragonal. In the first plane, O = /lua, and in 
the second O = —ha-^. Note that this means that for 
stripe order, there will be four ordering vectors, a pair at 
Q = ±\Q\{l + hg/ Ki^jO) from the first plane and a pair at 
Q = ±\Q\{0,1 + hg/ kl) from the other. However, for in- 
commensurate checkerboard order, there should be eight 
ordering vectors: ±\Q\{l + hg/KL,0), ±\Cl\{l — hg/KL,0), 
±|Q|(0, 1 + hg/KL), and ±|Q|(0, 1 - hg/nL). 



VI. ANALYSIS OF EXPERIMENTS IN THE 
CUPRATES 

There have been an extremely large number of ex- 
periments which have been performed on various closely 
related cuprates, both superconducting and not, which 
have been interpreted as evidence for or against the pres- 
ence of charge order of various types. For instance, 
there is a large amount of quasi-periodic structure ob- 
served in the local density of states measured by scan- 
ning tunneling microscopy (STM) on the surface of su- 
perconducting Bi2Sr2CaCu208+5 crystals, but there is 
controversy concerning how much of this structure arises 
from the interference patterns of well-defined quasiparti- 
cles whose dispersion is determined by the d-wave struc- 
ture of the superconducting gapSi&iSiiSi and how much 
refiects the presence of charge order or incipient charge 
order pi*2*^*^i2i& A similar debate has been carried out con- 
cerning the interpretation of the structures seen in inelas- 
tic neutron-scattering experimentsi^i^^i^°i^^i^^i^^i^'^i^^i^^ 

As mentioned in the introduction, the issue of how to 
distinguish charge order from interference patterns was 
discussed in detail in a recent review ji and so will not be 
analyzed here. Here, we will accept as a working hypoth- 
esis the notion that various observed structures should be 
interpreted in terms of actual or incipient order, and fo- 
cus on identifying the type of order involved. 



A. Neutron and X-ray scattering 

Scattering experiments in several of the cuprates, 
most notably La2-2;Sra;Cu04, Lai.6-z;Ndo.4Sra;Cu04, 
La2-2;Ba2:Cu04, and 0-doped La2Cu04 have produced 



clear and unambiguous evidence of charge and spin or- 
dering phenomena with a characteristic ordering vector 

which changes with doping|l°i27,28,29,30,31,32,33,34^35 ^he 

evidence is new peaks in the static structure factor corre- 
sponding to a spontaneous breaking of translational sym- 
metry, leading to a new periodicity longer than the lattice 
constant of the host crystal. In many cases, the period is 
near 4 lattice constants for the charge modulations and 
8 lattice constants for the spin. The peak-widths corre- 
spond to a correlation Iength2i2& that is often in excess 
of 20 periods. For technical reasons, the spin-peaks are 
easier to detect experimentally, but where both are seen, 
the charge ordering peaks are always seeni2i2L22i224£ to 
be aligned with the spin-ordering peaks, and the charge 
period is 1/2 the spin period. 

Except in the case^- of a very lightly doped (x < 
0.05) LSCO (where the stripes lie along an orthorhombic 
symmetry axis, so only two peaks are seen), there are four 
equivalent spin-ordering peaks and, where they have been 
detected, four equivalent charge ordering peaks. Thus, 
the issue arises whether this should be interpreted as the 
four peaks arising from some form of checkerboard order, 
or as two pairs of peaks arising from distinct domains of 
stripes - half the domains with the stripes oriented in the 
X direction and half where they are oriented along the 
y direction. A second issue that arises is whether the 
charge order is locked in to the commensurate period, 4, 
or whether it is incommensurate. 

A variety of arguments that the scattering pattern is 
revealing stripe order, and not checkerboard order, were 
presented in the original paper by Tranquada et. alm^ 
(and additionally in Ref. '38'44') where the existence of 
charge order in a cuprate high temperature supercon- 
ductor was first identified. Here, we list a few additional 
arguments based on the symmetry analysis performed in 
the present paper, which support this initial identifica- 
tion: 1. It follows from simple Landau theory^ that 
if there is non-spiral spin-order at wave- vectors Qi,Qj, 
there will necessarily be charge order at wave-vectors 
Qi + Qj- Thus, if the four spin-ordering peaks come from 
checkerboard order, then charge-ordering peaks should 
be seen at wave vectors ±2(5i, ±2(32 and ±Qi ± Q2, 
while if they come from stripe domains of the two ori- 
entations, no peaks at iQi ± Q2 should be seen. The 
latter situation applies to all cases in which charge order- 
ing peaks have been seen at all. 2. As mentioned above, 
in the LTT phase, the crystal fields should cause small 
splittings of the ordering vectors in an incommensurate 
checkerboard phase, causing there to be eight essentially 
equivalent Bragg peaks, as opposed to the four expected 
for domains of stripes of the two orientations. No such 
splittings have been detected in any of the scattering ex- 
periments on Lai.6-2;Ndo.4Sr2:Cu04 and La2_2;Ba2:Cu04 
crystals consistent with stripe domains. 3. It should 
be mentioned that the fact that the LTT phase stabi- 
lizes the charge order is, by itself, a strong piece of ev- 
idence that the underlying charge order is striped. In 
this phase, the O octahedra are tipped in orthogonal di- 



rections in alternating planes, and the direction of the 
tip is along the Cu-0 bond direction. This permits a 
uniquely strong coupling between the octahedral rota- 
tion and stripe order i^LSMSiii 

A second issue, especially when the period of the 
charge order is near 4 lattice constants, is whether the 
charge order is commensurate or incommensurate. One 
way to determine this is from the position of the Bragg 
peak - in the commensurate case, the structure factor 
should be peaked at 27r/4a (27r/8a for the spin order), 
and should be locked there, independent of temperature, 
pressure, or even doping for a finite range of doping. Most 
of the reported peaks seen in scattering are not quite 
equal to the commensurate value, however. In the LTT 
phase of La2_xBaa;Cu04, it is believed the stripe phase is 
locally commensurate. The ordering wave vector is tem- 
perature independent in the LTT phase, but jumps at the 
LTT-LTO transition and continues to change on warm- 
ing. For LSCO in the LTO phase, the stripes might be in- 
commensurate, however, there are only 4 peaks seen and 
not 8. So it must be incommensurate stripe order and not 
checkerboard order, A clearer piece of evidence comes 
from the rotation of the ordering vector away from the 
Cu-0 bond direction in the LTO phase of La2-a:Sr3:Cu04 
and O doped La2Cu04. In both cases, there is a small 
angle rotation (less than 4°) seen, which moreover de- 
creases with doping as the magnitude of the orthorhom- 
bic distortion decreases4& As discussed above, this is the 
generic behavior expected of incommensurate order, and 
is incompatible with commensurate order. 



B. STM 

The strongest quasiperiodic modulations seen in STM 
are those reported by Hanaguri et. alS^ on the surface of 
NaCCOC, which have a period which appears to be com- 
mensurate, 4 a. This observation has been interpreted as 
evidence that NaCCOC is charge-ordered with a checker- 
board pattern (at least at the surface^S) However, the 
correlation length deduced for the checkerboard order is 
only about two periods of the order. Indeed, the domain 
structure in the STM data looks to the eye very much 
like the pictures in our Figs. ^ (right panel) and |21 (right 
panel). This suggests the possibility that: 1: What is 
being seen is pinning of what, in the disorder free sys- 
tem, would be fluctuating order (a > 0) relatively close 
to the quantum critical point. 2: That the nearby or- 
dered state could be either a striped or a checkerboard 
state. We hope, in the near future, to apply the more 
quantitative analysis proposed in the present paper to 
this data. 

Concerning the modulations seen in STM studies on 
BSCCO: Given the recent interest in Bi2Sr2CaCu208+5, 
we report a preliminary application of our analysis to 
data from a near optimally doped sample, with an im- 
age size 21 CDW wavelengths across. Fig.[7|is a map of 
the LDOS integrated in energy to -|-15meVp^ (The axes 




FIG. 7: (color online) LDOS integrated in energy up to _B = 
+ 15meV. [Color scale is arbitrary.] Both S^orient and ^cdw 
are quite small, suggestive that the system is in a disorder- 
pinned, fluctuating phase. 



here are rotated 45° relative to those in Figs. and |21) 
In producing Fig. Elwe employ a Fourier mask (such as 
the one used Ref. \5w as a visual aid to show that there 
are indeed period 4 oscillations. This is a coherent state 
filter, centered in Fourier space around 27r/a(±l/4, 0) 
and 27r/a(0, ±1/4), and with a wide, flat top. Using the 
Fans. 111)14.71 we find ^orient = 4.5A and £,cdw = 2.5A, 
with A w 4.2a, and riorient — 0.28, which corresponds to 
7^1/2 and relatively strong disorder {Uq ~ 0.5). Ad- 
ditional measurements of the (unintegrated) LDOS on 
the same sample at E = SmeV, IbmeV yield compara- 
ble correlation lengths. From these we conclude the sys- 
tem shows a short-ranged mixture of (disorder-pinned) 
stripe and checkerboard order, and in the absence of 
pinning, would be in its fluctuating (symmetric) phase, 
but close to the critical point (a small). (Though there 
should probably be a fair amount of quasiparticle scat- 
tering at a nearby wave vector, it should be four-fold 
symmetric, so should not affect either ^orient or rjorient-) 
The fact that the orientational correlation length exceeds 
the CDW correlation length is suggestive that the prox- 
imate ordered state is a stripe ordered state and the ra- 
tio ^orient / DW ~ 2 is interesting, as it exceeds our 
(disorder-averaged) result of 1/2 for the symmetric phase 
(a > 0). However, undue weight should not be given to 
this result, as the (a > 0) region of Fig. [S] is a product 
of disorder-averaging, and Fig. [7| is a single set of data. 
In the future, we hope to apply our methods to a more 
substantial set of experimental data. 
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VII. CONCLUSIONS 

There are many circumstances in which charge order 
plays a significant role in the physics of electronically in- 
teresting materials. Depending on the situation, different 
aspects of the physics may be responsible for the choice 
of the characteristic period of the charge order; for in- 
stance, it can be determined by Fermi surface nesting (as 
in a Peierls transition), by a small deviation from a com- 
mensurate electron density (which fixes a concentration 
of discommensurations), or by some form of Coulomb 
frustrated phase separation. Working backwards, mea- 
surements of the period of the charge order as a function 
of parameters (temperature, pressure, doping, . . . ) can 
shed light on the mechanism of charge ordering. 

The physics that determines the ultimate pattern of 
charge order is still more subtle. For instance, for ad- 
sorbates on graphite, the sign of the energy of intersec- 
tion determines whether the discommensurations form a 
striped or honeycomb arrangement In 2H-TaSe2, bro- 
ken hexagonal symmetry has been observed-^'^ in x-ray 
scattering and TEM^'* (such a system has been studied 
by McMilliani^ using LG methods.) In certain nearly 
tetragonal rare-earth tellurides, which have been found 
to form stripe ordered phaseSf^^*^ this can be shown to 
be a consequence of some fairly general features of the 
geometry of the nested portions of the Fermi surface so 
long as the transition temperature is sufficiently highiS 

In the cuprates, calculations of the structures originat- 
ing from Coulomb- frustrated phase separatioupS, DMRG 
calculations on t-J laddersjS and Hartree-Fock calcu- 
lations on the Hubbard modelffiiSLS^ all suggest that 
stripe order is typically preferred over checkerboard or- 
der. Conversely, the Coulomb repulsion between dilute 
doped holes, or between dilute Cooper pairs favor a more 
isotropic (Wigner crystalline) arrangement of charges 
with more of a checkerboard structureiSiSMSiS& Thus, re- 
solving the nature of the preferred structure of the charge 
ordered states in the cuprates, at the least, teaches us 
something about the mechanism of charge ordering in 
these materials. 



On the basis of our present analysis, we feel that there 
is compelling evidence that most, and possibly all, of 
the charge order and incipient charge order seen in hole- 
doped cuprates is preferentially striped. We also con- 
clude that most of the structure seen in STM studies is 
disorder pinned versions of what would, in the clean limit, 
be fiuctuating stripes, rather than true, static stripe or- 
der. 

Note: After this work was completed we received a 
draft of a paper by del Maestro and coworkers^ who 
discuss similar ideas to the ones we present in this paper. 
We thank these authors for sharing their work with us 
prior to publication. After this paper was submitted for 
publication M. Vojta pointed out to us that in a very 
recent paper he and his coworkers considered the effects 
of slow thermal fluctuations of stripe and checkerboard 
charge order on the magnetic susceptibility of disorder- 
free high Tc cuprates. "S4. 

Note added: While this paper was being refereed a new 
neutron scattering study of LNSCO became availableSS, 
which confirmed the existence of unidirectional charge 
order (stripe) and coUinear spin order in this material, in 
agreement with the results and interpretation of Ref. [l^ . 
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